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Abstract 


A simple method is presented to estimate the complex dielectric constants of 
individual layers of a multilayer composite material. The multilayer composite material 
sample is loaded in an X-band rectangular waveguide and its two port 5-parameters are 
measured as a function of frequency using the Hewlett-Packard (HP) 8510 Network 
Analyzer. Also, by applying the mode matching technique, expressions for the S- 
parameters of the composite material as a function of electric properties of individual 
layers are developed. Using the MatLab Optimization Tools simple MatLab scripts are 
written to search for electric properties of individual layers so as to match the measured 
and calculated S-parameters. 

A single layer composite material formed by using materials such as Bakelite, 
Nomex Felt, Fiber Glass, Woven Composite B and G, Nano Material #0, Cork, Garlock, 
of different thicknesses are tested using the present approach. The dielectric constants of 
these materials estimated using the present approach (assuming the thicknesses are 
known) are in good agreement with their true values. Assuming the thicknesses of 
samples unknown, the present approach is shown to work well in estimating the 
dielectric constants and the thicknesses. 

A number of two layer composite materials formed by various combinations of 
above individual materials are tested using the present approach. The values of dielectric 
constants of individual layers estimated (assuming the thickness of individual layers are 
known) using the present approach are in good agreement with their true values. 

However, the present approach could not provide estimate values close to their true 
values when the thicknesses of individual layers were assumed to be unknown. This is 
attributed to the difficulty in modelling the presence of airgaps between the layers while 
doing the measurement of S-parameters. A few example of three layer composite are 
also presented. 


I. Introduction 

Multilayer substrates are used for many practical applications such as 
Microwave Integrated Circuits (MIC), Monolithic MIC ( MMIC) [1] , radomes for 
protection of antennas from external environment, spatial filters for antenna beam 
shaping [2], and Frequency Selective Surfaces (FSS) [3-4]. With proper choices of 
individual layers in a multilayer composite substrates it is possible to achieve a 
composite material with altogether new properties that were otherwise not found in the 
individual layers. Exact knowledge of the material properties such as permittivity and 
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permeability of individual layers in a multilayer substrate is essential for designing spatial 
filters, radomes and composite materials for FSS applications. However, the present 
practice in the estimation of electric properties of composite materials with multi-layers 
emphasizes only the determination of overall effective properties. In this report an 
attempt is made to provide a procedure to estimate the electric properties of individual 
layers of a composite material. 

Permittivity and permeability of these composite multi-layer substrates can be 
estimated using one of the following methods: 1) free-space techniques; 2) resonant 
cavity perturbation techniques; and 3) transmission line methods. Each method has its 
own advantages and limitations. For measurement of electric properties over a wide 
frequency range, transmission line or waveguide methods [5-8] are more popular even 
though they are less accurate due to unavoidable measurement errors. In the waveguide 
measurement method, a sample of composite multilayer substrate is placed in a 
waveguide and two port S-parameters are measured using a hp-8510 Network Analyzer. 

In the earlier work [7-8], assuming that an equivalent homogeneous material occupies 
the sample space, 5-parameters are estimated as a function of effective electric properties 
of the equivalent material. Using the inverse procedure, constituent parameters of the 
equivalent material are determined by matching the estimated and measured 5- 
parameters. In these methods electric properties of individual layers are not determined. 
However, for some applications such as radomes and FSS, a knowledge of electric 
properties of individual layers is required. In this report, a waveguide measurement 
method is presented to estimate the electric properties of individual layers of a composite 
multilayer substrate. 
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The remainder of this report is organized as follows. The mode matching 
formulation of a waveguide loaded with multilayer composite material is developed in 
section H. Also in section n, the expression are developed to estimate 2-port S- 
parameters of composite material slab. Numerical and measured results on the S- 
parameters of numerous composite material are presented in section HI for the direct 
problem where electrical properties of individual layers are assumed to be known. Also 
in section HI, using the MatLab Codes given in Appendices, estimated values of 
dielectric constants of one, two, and three layer composite material are presented. The 
report concludes in section IV with remarks on the validity and usefulness of the present 
method. 

II. Theory 

In this section, the method of moments is used to determine the S-parameters of 
a rectangular waveguide loaded with a multilayer dielectric substrate as shown in figure 
1 . 


Y Incident Wave p or t #\ 



Figure 1 : Geometry of rectangular waveguide loaded with 
multilayer composite material. 


3 



The multilayer substrate consists of N layers, where the n' h layer has parameters 
) and is located between the transverse planes at z — z„_, and z = z„. 

S,, and S 21 Parameters: 

To estimate the S u and S 21 components of 2-port parameters, it assumed that the 
TE l0 mode of unit amplitude is incident on the interface at z = 0 from the region z <0 . 

If E w ,E n ....E tN are the transverse electric fields on the interfaces at z = 0 ,z l ,Z 2 '<—Z N , 

respectively, then the transverse electric fields in the various regions of waveguide are 
obtained as 
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The transverse electric fields over the interfaces can be expressed in terms of vector 
modal expansion functions as: 
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Equations (1 1) -(13) are the required integral equations to be used to determine the 
complex amplitudes T 0j ,T tj ,....T Nj . Selecting h k as a testing function and using the 

Galerkin’s procedure, the equations (1 1)-(13) are converted into a set of simultaneous 
equations: 
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Due to the orthogonal nature of vector modal functions it can be shown that the complex 
amplitudes T 0J ,j = 1,2,3. ..7 0 , T w j = 1,2,3..../, , 7^, j = 1,2,3... J N are all zeros. 

Hence equations (14)-(16) can be simplified as 
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solution of above (N + 1) equations gives an estimate of complex amplitudes 
^oo’^io’Ta)’ T n0 from which 5,, and S 21 are determined as 


S„ =T 0 o-l and S 2l = T m e iPl ' N 


(18) 


S 22 and 5 12 Parameters: 

The port 2 parameters, S 22 and 5 I2 can be determined by following the procedure used 
for estimation of 5 U and 5 21 , and reversing the locations of the layers as shown in 
figure 2. Note that the S 22 calculated using the reference planes shown in figure 2 and 

S 22 measured using the hp-8510 network analyzer differ by phase e l/s ° 1ZN . 



Figure 2: Geometry of rectangular waveguide loaded 
with a composite material for estimation of S 22 and 5 12 
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III. Numerical Results 


A: Estimation of S-Parameters( Direct/Forward Problem) 


A simple MatLab code (Appendix A) is written to solve the simultaneous equations 
given in equations (17) and determine all four S -parameters of a composite material slab 
placed in a rectangular waveguide. In this section, assuming the properties of individual 
layers of a composite material known, the S -parameters for various composite slabs are 
computed (using the MatLab code) as a function of frequency and compared with the 
measured S -parameters. 

Single Layer Composite Material: 

Figure 3 shows S u and S 2l parameters of a composite material consisting of a 
Garlock single layer as a function of frequency. In this case a single layer of Garlock 
with thickness A, = 0. 17cm and electric properties £ r = 7.5 - y'0.001 is used to form a 
composite material. Excellent agreement between measured and estimated 


tthgla QirlBcfciajroftlTlBlfi— a 17 cm 

— 



Figure 3: Measured and estimated S -parameters of single 
Garlock slab. Thickness Aj = 0.17cm, £ r =1.5- y‘0.001, 
ju r =l.0-j0.0 
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values of the S -parameters validates the MatLab code. Figure 4 shows the S u and S 2I 
parameters of a composite material formed by a single layer of Teflon material 
( e r =2.03- y'0.001, jd r =1.0- jO.O ) of thickness A, =0.635 cm. A good agreement 
between the measured and estimated values of S -parameters confirms validity of the 
present method. Note that the parameters S 22 and S l2 are expected and found to be 
identical to 5,, and S 21 , respectively. 


HiieniTtotWd w aiMgii 

i n R«||(B11) 

— — l mag i ) ca*umm 



Figure 4: Measured and estimated S -parameters of single teflon slab. 
Thickness A, = 0.635cm , £ r =2.03 — j‘0.001, ju r =1.0- jO.O 


Two Layer Composite Material: 


For further validation of the present method and MatLab code, the S -parameters of 


composite material formed by various combination of two layers are fabricated and 


tested. A first sample considered consists of Bakelite and Teflon layers. 
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The composite slab is formed by placing the bakelite layer with A, = 0.33cm , 


z 0 =0, z, =0.33 cm and the Teflon layer with A 2 = 0.635cm, z, = 0.33cm, 
z 2 = 0.965cm . Measured and estimated S u and S 21 parameters as a function of 
frequency are shown in figure 5. For the composite material described in figure 5, the 
measured and estimated S 22 and S l2 values are shown in figure 6. 


BakeHte -Teflon layers 

Bake lie d = 0.3302 cm,er = 3.76 -1 0.001 

Teflon d = 0.635 cm, er = 2.03 -J0.Q01 
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Figure 5: Measured and estimated S -parameters of two layers composite 
material (Bakelite-Teflon). Bakelite: A, = 0.33cm, z 0 =0, z, = 0.33cm, 

£ r =3.76- j0.00l,ju r =1.0, Teflon: A, = 0.635cm, z, = 0.33cm, 

Z 2 =0.965 cm, £ r = 2.03- j0.00l,ju r =1 
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Bate Ito- Tftflon laymr a 

Bate He d = 0.3302 cm, or = 3.76 -JO .001 

Tftflon d = 0.635 cm ar = 2.03-J0.001 


Real(S22) 

- — — - lmag.(S22) 



Figure 6: Measured and estimated values of S 22 and S l2 for composite 
slab (parameters as described in figure 5) 


A second sample of composite material considered consists of Garlock and Nomex 
materials. The composite slab is formed by placing first a Garlock slab of thickess 
A, = 0.17cm, e r = 7.5- y'0.001,// r =1.0 at z 0 =0 and z, = 0.17cm and then a Nomex 

slab of thickness A 2 = 0.33cm , e r = 1.2 - j0.00l,ju r =1.0 at z, = 0.17cm and 

Z 2 = 0.5cm . Measured and estimated A,, , 5 2 , , 5 12 , and S 22 parameters for Garlock- 

Nomex composite slab are shown in figures 7 and 8. 

From figures 7-8, estimated values of S -parameters for the Garlock-Nomex composite 
slab agrees well with measured values of the S -parameters. However, for the Bakelite- 
Teflon combination there is small disagreement between measured and estimated S 22 and 
S',, parameters. This disagreement may be attributed to the presence of an air gap 
between the slabs of Bakelite and Teflon. 
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Figure 7: Measured and estimated S -parameters of two layers 
composite material (Garlock-Nomex). Galock: =0.17 cm, z 0 =0, 

Z, = 0.17cm, e r = 7.5- j'0.001,// r =1.0,Nomex A 2 =0.33 cm, 

Z| = 0.17cm, z, = 0.50cm, = 1 . 2 - 70 . 001 ,//,. =1 



Figure 8: Measured and estimated S-parameters of two layers composite 
material (Garlock-Nomex) with dimensions as shown in Figure 7 
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Three Layer Composite Material: 


A third sample of composite material considered consists of Nano Material # 0, 
Garlock and Garlock slabs. The composite slab is formed by placing first a Nano 
Material #0 slab of thickess A l = 0.31cm , £ r = 2.5 - j0.00l,jU r =1.0 at z 0 = 0 and 

Z, = 0.31cm and then a Garlock slab of thickness A 2 = 0.17cm , 

£ r =7.5- j0.00l,ju r =1, at Zi = -31cm, and z 2 = 0.48cm. Another Garlock slab of 
thickness A, = 0.17cm, £ r =7.5-y'0.001,// r =1.0 is placed at z 2 = 0.48cm and 
z 2 = 0.65cm . Figures 9-10 show measured and estimated S -parameters of composite 
slab consisting of Nano_Mat #0-Garlock-Garlock combination. 


Composite Slab: Nano Mateiil #0, Garlock, Garlock 

Real(SII) 



Figure 9: Measured and estimated 5,, and S 2l parameters of composite slab 
consisting of Nano-Material#0, Garlock, and Garlock(Dimensions as 
described in text) 
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Composite Slab: Nano Materil #0, Garlock, Garlock 



Figure 10: Measured and estimated S 22 and S 12 parameters of 
composite slab consisting of Nano-Material#0, Garlock, and 
Garlock(Dimensions as described in text) 


From figures 9-10, the trend in the values of S -parameters predicted by the computation 
is similar to the trend observed in the measured values. However, significant differences 
are observed between the measured and estimated values of S -parameters. This may be 
attributed to the sample preparation where the presence of air gaps between the 
individual slabs is unavoidable. Also even though these individual slabs may be in 
physical contact with each other, electrical contact may still not be insured. It is also 
noticeable that the disagreement between the measured and estimated S -parameters gets 
worse as more individual slabs are used to form a composite slab. 
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B: Estimation of Dielectric Constants( Inverse Problem): 

Single Layer Composite Material: 

In this section a procedure is described for the computation of the complex 
dielectric constant of a given composite sample from the two port measured data. For a 
given composite material sample, let the two port S -parameters be (S llm , S 2im , 

S| 2 m , S 22m ), measured using the hp-8510 Network Analyzer. With a prior knowledge of 
the thickness of composite slabs, the two port S -parameters can be estimated using the 
MatLab Code given in Appendix A as a function of £ r . Let the estimated values of 5 - 
parameters be 5, lc (£ r ), S 2lc (£ r ),S l2c (£ r ), S 22c (£ r ) . The errors in estimated and measured 
S -parameters can then be written as ER { = real(S Uc - S llm ) , ER 2 = imag.(S llc — S„ m ) , 

ER 2 = real(S 2lc - S 2lm ) , ER 4 = imag.(S 2lc - S 2lm ) , ER n = real(S 22c - S 22m ) , 

ER S = imag.(S 22c - S 22m ) . The total mean squared error or the objective function to be 
minimized as a function of £ r can be written as 

ER t = -yJiER; + ER; + ER; + ER; + ER; + ER] + ER: j + ER ; ) (19) 

A simple MatLab Code given in Appendix B minimizes the objective function in (19) to 
estimate the unknown value of dielectric constant £ r . 

Examples (When Thickness of Sample Is Known): 

Garlock Slab: A single layer of Garlock slab of size (2.29 x 1.02 x 0.17 )cm was placed 
in an X-band rectangular waveguide. After proper calibration of the hp 8510 Network 
Analyzer, two port S-parameters of the Garlock slab were measured as a function of 
frequency. The measured data is stored in a file FT_GRLK.60. Using the MatLab Code 
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given in the Appendix B, the dielectric constant of the Garlock slab is estimated and 


shown in figure 1 1 . 


Real and Imaginary Parts of Relative Dielectric 
Constant of Garlock Material 
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Figure 1 1 : Relative dielectric constant (real and imaginary) of 
Garlock material estimated using measured and ideal S- 
parameters. 


In figure 1 1 solid and dash-dot lines are the estimates of dielectric constant using the 
ideal or noise free S-parameters calculated assuming the dielectric constant is known and 
equal to e r = 7.5 - / 0.001 . The estimates shown in figure 1 1 were obtained by 
considering the thickness of the slab equal to 0.170cm. To check the level of confidence 
in these estimates, from the estimated values of dielectric constant (estimated using the 
measured data), the difference between measured and computed values of S-parameters 
is plotted in Figure 12. 
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Figure 12: Differences between computed and measured values of S- 
parameters. Calculated S-parameters are determined using estimated 
dielectric constant 


The differences between the computed and measured S-parameters are within the limits 
set in the optimizer. Ideally, the difference between the computed and measured S- 
parameters must be close to zero. However, to achieve ideal results, the optimizer would 
take longer time. 

Number of single layer composite materials of various material and thickness were 
constructed and their S-parameters were measured over the X-band frequency range. 
From these measured values of S-parameters and using the MatLab Code given in 
Appendix B, the dielectric constants of these single layer composite were estimated and 
presented in figures 13-25. 
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Real and Imaginary Parts ot Relative Dielectric 
Constant ol Nfano-Mat 0 Mater tat 



Figure 13: Relative dielectric constant (real and imaginary) of Nano 
material #0 estimated using measured and ideal S-parameters 
(thickness = 0.3099 cm) 


Real and Imaginary Parts of Relative Dielectric 
Constant of Teflon Material 



Figure 14: Relative dielectric constant (real and imaginary) of Teflon 
material estimated using measured and ideal S-parameters (thickness 
0.9398 ) 
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Real and Imaginary Parts of Relative Dielectric 
Constant of Cork Material 



Frequency(GHz) 


Figure 15: Relative dielectric constant (real and imaginary) of Cork 
material estimated using measured and ideal S-parameters (thickness 
0.3048cm) 


Real and Imaginary Parts of Relative Dielectric 
Constant of Ceramic Material 



Figure 16: Relative dielectric constant (real and imaginary) of Ceramic 
material estimated using measured and ideal S-parameters (thickness = 
0.2845 cm) 
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Real and Imaging Parts of Relative Dielectric 
Constant of Nomex Felt Material 



Figure 17: Relative dielectric constant (real and imaginary) of 
Nomex Felt material estimated using measured and ideal S- 
parameters (thickness = 0.33cm) 
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Figure 18: Relative dielectric constant (real and imaginary) of Rubylith 
material estimated using measured and ideal S-parameters (thickness = 
0.3683 cm) 
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Real and Imaginary Parts of Relative Dielectric 
Constant of Bakelite Material 



Frequency(GHz) 

Figure 19: Relative dielectric constant (real and imaginary) of Bakelite 
material estimated using measured and ideal S-parameters (thickness = 
0.33cm 


Real and Imaginary Parts of Relative Dielectric 
Constant of Fiber Glass Material 



Frequency(GHz) 


Figure 20: Relative dielectric constant (real and imaginary) of Fiber 
Glass material estimated using measured and ideal S-parameters 
(thickness = 0.0533cm 
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Real and Imaginary Parts ol Relative Dielectric 
Constant of Woven Composite B Matenal 



Figure 21: Relative dielectric constant (real and imaginary) of Woven 
Composite B material estimated using measured and ideal S-parameters 
(thickness = 0.1625cm 


Real and Imaginary Parts of Relative Dielectric 
Constant of Low Density Foam Material 



Figure 22: Relative dielectric constant (real and imaginary) of low density Foam 
material estimated using measured and ideal S-parameters (thickness = 
0.2921cm) 
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Real and imaginaty Parts of Relative Dielectric 
Constant of Woven Composite G Material 



Figure 23: Relative dielectric constant (real and imaginary) of Woven 
Composite G material estimated using measured and ideal S-parameters 
(thickness = 0.2083 cm) 

Examples (When Thickness of Sample Is Unknown): 

In the previous section it was assumed that the thickness of the sample is known 
apriori. These thicknesses are measured in the laboratory using an electronic 
micrometer. However, for compressible and thin samples the accuracy of these 
measurements is questionable. In this section it is assumed that the thickness of the 
sample is unknown and the optimizer is asked to estimate the thickness along with the 
dielectric constant of a single layer composite material slab. The MatLab Code used for 
estimation of sample thickness as well as dielectric constant is given in Appendix B. 
Garlock Slab: The measured data stored in a file FT_GRLK.60 is used to estimate the 
thickness of Garlock slab and dielectric constants. Using the MatLab Code given in the 
Appendix B, the dielectric constant of the Garlock slab and its thickness are estimated 
and shown in figure 24. 
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Real and Imaginary Parts of Relative Dielectric 
Constaht of Garlock Material (Thickness Unknown) 



Figure 24: Relative dielectric constant (real and imaginary) of Garlock 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.17 cm. Thin 
solid line indicate estimated value of sample thickness. 


In figure 24, estimated values of the dielectric constant assuming sample thickness 
unknown are in close agreement with the dielectric constant estimated assuming known 
value of sample thickness. The average value of estimated sample thickness is 0.1699 cm 
which is very close to the actual measured thickness of 0.1702cm. This validates the idea 
that apriori knowledge of sample thickness is not necessary for the inverse problem. In 
fact the thickness of the sample can be considered as one of the unknown variables along 
with the dielectric constants to optimize the error function defined in (19). Using the 
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Figure 25: Differences between computed and measured values of S- 
parameters. Calculated S-parameters are determined using estimated 
dielectric constant and thickness of the slab 


estimated values of dielectric constants and the sample thickness, the difference between 
measured and computed values of S-parameters is plotted in Figure 25. 

The differences between the computed and measured S-parameters are within the limits 
set in the optimizer. Ideally, the difference between the computed and measured S- 
parameters must be close to zero. However, to achieve ideal results, the optimizer would 
take a longer time. 

From the measured values of the S-parameters for a variety of single layer 
composite materials (Nano Material #0, Cork, Ceramic, Nomex Felt, Rubylith, Bakelite, 
Teflon, Fiber Glass, Woven Composite B and G, Nano Material # 0) the dielectric 
constants and thicknesses of the material are estimated using the MatLab Code given in 
Appendix B. These estimated values are shown in figures 26-36. For comparison, the 
estimated values of dielectric constants using MatLab Code given in Appendix B are also 
plotted in figures 26-36. The estimated values of dielectric constants assuming the 
thickness unknown are in good agreement with the estimates obtained using known 
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values of dielectric constants. From the figures 26-36 it may be concluded that the 
thickness of a single layer slab can be treated as one of the unknown variables along with 
the dielectric constants. 


Real and Imaginary Parts of Relative Dielectnc 
Constant of Nano Material #0 (Thickness Unknown) 
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Figure 26: Relative dielectric constant (real and imaginary) of Nano 
material #0 estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.301 cm. Thin 
solid line indicate estimated value of sample thickness. 
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Real and Imaginary Parts of Relative Dielectric 
Constant of Cork Material (Thickness Unknown) 
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Figure 27: Relative dielectric constant (real and imaginary) of Cork 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.3048 cm. 
Thin solid line indicate estimated value of sample thickness. 


Real and Imaginary Parts of Relative Dielectric 
Constant of Ceramic Material (Thickness Unknown) 



Figure 28: Relative dielectric constant (real and imaginary) of Ceramic 
material estimated assuming thickness unknown. Symbols indicate estimated 
values of dielectric constant using thickness = 0.2845 cm. Thin solid line 
indicate estimated value of sample thickness. 
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Real and Imaginary Parts of Relative Dielectric 
Constant of Nlomex Felt Material (Thickness Unknown) 
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Figure 29: Relative dielectric constant (real and imaginary) of Nomex Felt 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.3302 cm. Thin 
solid line indicate estimated value of sample thickness. 


Rea) and Imaginary Parts of Relative Dielectric 
Constant ol Rubylith Material (Thickness Unknown) 
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Figure 30: Relative dielectric constant (real and imaginary) of Rubylith 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.3683 cm. 
Thin solid line indicate estimated value of sample thickness. 
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Real and Imaginary Parts of Relative Dielectric 
Constant of E&kelrte Material (Thickness Unknown) 
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Figure 3 1 : Relative dielectric constant (real and imaginary) of Bakelite 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.3302 cm. Thin 
solid line indicate estimated value of sample thickness. 


Real and Imaginary Parts of Relative Dielectric 
Constant of Teflon Material (Thickness Unknown) 
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Figure 32: Relative dielectric constant (real and imaginary) of 
Teflon material estimated assuming thickness unknown. Symbols 
indicate estimated values of dielectric constant using thickness = 
0.9398 cm. Thin solid line indicate estimated value of sample 
thickness. 
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Real and Imaginary Parts of Relative Dielectric 
Constant of Fiber Glass Material (Thickness Unknown) 



Figure 33: Relative dielectric constant (real and imaginary) of Fiber Glass 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.0533 cm. Thin 
solid line indicate estimated value of sample thickness. 
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Figure 34: Relative dielectric constant (real and imaginary) of Foam 
material estimated assuming thickness unknown. Symbols indicate 
estimated values of dielectric constant using thickness = 0.2921 cm. 
Thin solid line indicate estimated value of sample thickness. 
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Real and Imaginary Parts of Relative Dielectric 

Constant of Woven Composite B Material (Thickness Unknown) 
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Figure 35: Relative dielectric constant (real and imaginary) of Woven Composite 
B material estimated assuming thickness unknown. Symbols indicate estimated 
values of dielectric constant using thickness = 0.1626 cm. Thin solid line 
indicate estimated value of sample thickness. 


Real and Imaginary Parts of Relative Dielectric 

Constant of Woven Composite G Material (Thickness Unknown) 



Frequency(GHz) 

Figure 36: Relative dielectric constant (real and imaginary) of Woven 
Composite G material estimated assuming thickness unknown. Symbols 
indicate estimated values of dielectric constant using thickness = 0.2083 cm. 
Thin solid line indicate estimated value of sample thickness. 
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Most of estimated dielectric constants determined assuming thickness unknown agree 


well with the dielectric constants estimated using measured value of slab thickness. 
However, in figure 33, the two results disagree significantly. This is because the 
thickness of the slab is very small. In fact, the S-parameters calculated using estimated 
thickness and the dielectric constants agree very well with the measurements as shown in 
figure 37. 


Skigla Fber Glass layer (TNcknftSt Unknown) 



Froquency(GHz) 

Figure 37: Measured and computed S-parameters for Fiber Glass Material. 
Computed S-parameters are determined using estimated values of dielectric 
constant and thickness of Fiber Glass. 


Two Layer Composite Material: 

In this section measured S-parameters of composite material consisting of two material 
layers are used to estimate the dielectric constants of individual layers. The error or 
objective function used for this purpose is identical to equation (19). The MatLab Code 
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given in Appendix B, with proper input variables is used to estimate the dielectric 
constants of two layers. 

Examples: (Thickness of Layers Known): 

A two layer composite material consisting of Bakelite and Teflon was formed by placing 
first the Bakelite layer of thickness A, = 0.33 cm between the z 0 = 0 and z, = 0.33 cm 

planes. A Teflon layer of thickness A 2 = 0.635 cm was then placed between z, = 0.33cm 
and z 2 = 0.965 cm . The measured values of S-parameters of the Bakelite-Teflon as a 
function of frequency was stored in a file FT_BKLTEF.60. Using the MatLab Code 
given in Appendix B, the dielectric constants of the two layers were estimated and 
plotted in Figure 38. 


Two Layer* Composite Material : Bakettte-Teflon 
(Thicknesses of layers assumed known) 



Figure 38: Relative dielectric constant (real and imaginary) of Bakelite-Teflon 
Composite material estimated using measured S-parameters (assuming thickness 
known) 
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Figure 39 shows estimated values of dielectric constants of Bakelite-Teflon Composite 
material using the ideal S-parameter values computed using the MatLab Code given in 
Appendix A. Figure 40 shows the error involved in the estimations of dielectric 
constants using measured and ideal values of S-parameters. From figure 40 it is clear that 
the error in estimation using the measured S-parameters is higher than the error involved 
in the estimation using ideal S-parameters. This may be attributed to the noise present in 
the measured data and inability of the analytical model to take into account potential air 
gaps present between the two layers. 
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Two Layers Composite Material : Bakelite-Teflon 
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Estimates are determined using ideal S-parameters 
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Figure 39: Relative dielectric constant (real and imaginary) of Bakelite-Teflon 
Composite material estimated using ideal S-parameters (assuming thickness known) 
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Optimum Value o t Error Function 



Figure 40: Value of error function given in (19) in estimating the dielectic 
constants of Bakelite-Teflon Composite using both measured and ideal S- 
parameters 


For further verifications of the present approach for a two layer composite material the 
following samples were considered: 

1) Garlock-Nomex Felt: A, = 0.17cm, A, = 0.33cm , z 0 =0.0, z, = 0. 1 1cm , 

Z 2 =0.5 cm 

2) Garlock-Ceramic: A, = 0.17cm, A, = 0.2845cm, z 0 =0.0, z, = 0.17cm, 

Z 2 = 0.4545cm 

3) Garlock-Nano Material #0 : A, =0.1 7cm, A, = 0.3099cm, z 0 =0.0, 
z,=0.17cm, z 2 = 0.4799cm 
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4) Garlock-Woven Composite G: A, = 0.17cm , A, = 0.2083cm , z 0 = 0.0, 

Z, = 0.17cm, z 2 - 0.3783cm 

5) Garlock-Bakelite A, = 0.17cm , A 2 = 0.33c/?; , z 0 =0.0, z, — 0.17cm, 

Z 2 = 0.50cm . 

For these samples the relative dielectric constants of individual layer were estimated 
from the measured and ideal S-parameters. Estimated values of relative dielectric 
constants are shown in figures 41-56. 
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two Layers compost* Material : (ianock-Nomex hett 
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Figure 41: Relative dielectric constant (real and imaginary) of Garlock-Nomex 
Felt Composite material estimated using measured S-parameters (assuming 
thickness known) 
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Two Layers Composite Material : Garlock- Nomex Felt 
(Thicknesses of layers assumed known) 

Estimates are determined using ideal S-parameters 
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Figure 42: Relative dielectric constant (real and imaginary) of Garlock - Nomex Felt 
Composite material estimated using ideal S-parameters (assuming thickness known) 


Optimum Value of Error Function 



Figure 43: Value of error function given in (19) in estimating the dielectic constants of 
Garlock -Nomex Felt Composite using both measured and ideal S-parameters. 
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Two Lay«r« Composite Material : Garlock-Coramic 
(Thlctatesass of layer* assumed known) 



Frequency (GHz) 

Figure 44: Relative dielectric constant (real and imaginary) of Garlock-Ceramic 
Composite material estimated using measured S-parameters (assuming thickness 
known) 
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Figure 45: Relative dielectric constant (real and imaginary) of Garlock - Ceramic 


Composite material estimated using ideal S-parameters (assuming thickness known) 
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Optimum Value of Error Function 



Frequency (GHz) 


Figure 46: Value of error function given in (19) in estimating the dielectic constants of 
Garlock -Ceramic Composite using both measured and ideal S-parameters. 


Two Layers Composite Material : Gartocfc-Nano Material #0 
(Thtcknaeeae of leyere assumed known) 



Figure 47: Relative dielectric constant (real and imaginary) of Garlock-Nano Material #0 
Composite material estimated using measured S- parameters (assuming thickness known) 
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Two Layers Composite Material : G a rlock -Nano Material #0 
(Thicknesses of layers assumed known) 

Estimates are determined using Ideal S parameters 
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Figure 48: Relative dielectric constant (real and imaginary) of Gariock - Nano 
Material #0 Composite material estimated using ideal S-parameters (assuming 
thickness known) 
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Figure 49: Value of error function given in (19) in estimating the dielectic constants 
of Gariock -Nano Material #0 Composite using both measured and ideal S-parameters. 


42 


Two Layers Composite Material : Garlock-Woven Composite G 
(Thicknesses of layers assumed known) 



Frequency (GHz) 

Figure 50: Relative dielectric constant (real and imaginary) of Garlock-Woven 
Composite G material estimated using measured S-parameters (assuming thickness 
known) 


Two Layers Composite Material : Garlock-Woven Composite G 
(Thicknesses of layers assumed known) 

Estimates are determined using Ideal S-parameters 


d> 

co 


o 

cc 

sr 

c 

a 


M 


S 


o 


a 

a 

a 



Frequency (GHz) 


Figure 51: Relative dielectric constant (real and imaginary) of Garlock - Woven 
Composite G material estimated using ideal S-parameters (assuming thickness known) 
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Optimum Value of Error Function 



Figure 52: Value of error function given in (19) in estimating the dielectic constants of 
Garlock -Woven Composite G material using both measured and ideal S-parameters. 


Two Layers Composite Material : Gartock-Bakatfte 
(Thicknesses of layers assumed known) 



Figure 53: Relative dielectric constant (real and imaginary) of Garlock-Bakelite 
Composite material estimated using measured S-parameters (assuming thickness 
known) 
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Two Layers Composite Material : Gartock-Bakofite 
(Thicknesses of layers assumed known) 

Estknates are determined using ideal S-parameters 
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Figure 54: Relative dielectric constant (real and imaginary) of Garlock - Bakelite 
Composite material estimated using ideal S-parameters (assuming thickness 
known) 

Opttnum Value at £ mar Function 



Figure 55: Value of error function given in (19) in estimating the dielectic 
constants of Garlock -Bakelite Composite material using both measured and 
ideal S-parameters. 
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From the results shown in figures 41-55, it may be concluded that the method presented 
in this report can successfully estimate dielectric constants of two layer composite 
material with apriori knowledge of the thickness of individual layers. An attempt was 
made to estimate thicknesses of individual layers along with their dielectric constants 
using the MatLab Code given in Appendix B. However, the optimizer did not estimate 
these values correctly and hence the results for such cases are not presented. 

Three Layer Composite Material: 

In this section measured S-parameters of composite material consisting of three material 
layers are used to estimate dielectric constants of individual layers. The error or 
objective function used for this purpose is identical to equation (19). The MatLab Code 
given in Appendix B, with proper input variables is used to estimate the dielectric 
constants of two layers. 

Examples:(Thickness of Layers Known): 

A three layer composite material consisting of layers of Nano Material #0, Garlock, and 

Garlock was formed by placing the Nano Material #0 layer of thickness 

A, = 0.3099cm between z 0 = 0 and z, = 0.3099cm planes, the Garlock layer of 

thickness A, = 0.17cm between z { - 0.3099cm and z 2 = 0.4799cm planes, and 
Garlock layer of thickness A, = 0.17cm between z 2 =0.4799 cm and z 3 = 0.6499cm 

planes. The measured values of S-parameters stored in FT_NANGRKGRK.60 were used 
to estimate dielectric constants of each individual layers using the MatLab code given in 
the Appendix B. The estimated values of dielectric constants are shown in Figure 56. 

The values of dielectric constants of each these material when estimated using single 
layer measured data were found to be £ r = 2.5 - j 0.00 for the Nano Material #0, 
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£ r = 7.5 - y'0.00 for the Garlock material. However, when the individual dielectric 
constants are estimated using the measured data for the three layer composite material, 
these estimates come out to be little different for the earlier estimates. This may be 
attributed to the presence of airgaps between the layers. 
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Figure 56: Estimated values of dielectric constants of individual layers from 
measured values of S-parameters of three layer composite slab (Nano Material 
#0-Garlock-Garlock) 



Figure 57: Final values of error function in estimation of dielectric constants of 
individual layers shown in Figure 56 
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The final value of error function or objective function for the estimation shown in Figure 
56 is shown in Figure 57. 

A second example of three layer composite material considered was formed by three 
layers of the same Garlock slabs. The first Garlock slab of thickness A, = 0.17cm was 
placed between the z 0 =0 and z, = 0.17cm planes, second Garlock layer of thickness 
A, = 0.17cm was placed between z { = 0.17cm and z 2 =0.34 cm planes, and the third 
Garlock layer of thickness A, = 0.17cm between z 2 = 0.34cm and z 3 =0.51 cm planes. 

The measured values of S-parameters stored in FT_GRKGRKGRK.60 were used to 
estimate dielectric constants of each individual layers using the MatLab code given in 
Appendix B. The estimated values of dielectric constants are shown in Figure 58. The 
values of estimated dielectric constants of each individual layers must be close to 
£ r = 7.5 - y'0.00 . However, the results shown in the Figure 58 are widely spread around 
the correct value of the dielectric constant. This discrepancy in the estimated values 
derived using the measured S-parameters of multilayer composite material is attributed to 
the measurement errors caused by the airgaps between the layers. The error function or 
objective function at the final values of dielectric estimates is shown in Figure 59. 
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Figure 58: Estimated values of dielectric constants of individual layers from 
measured values of S-parameters of three layer composite slab (Garlock-Garlock- 
Garlock) 



Figure 59: Final values of error function in estimation of dielectric constants 
of individual layers shown in Figure 58 
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IV Conclusion 

A simple waveguide mode matching method in conjunction with th tfminsearch 
MatLab optimization search function has been presented to estimate the complex 
permitivity of multi-layer composite material. A multi-layer composite material is placed 
in an X-band rectangular waveguide and its two-port S-parameters are measured over the 
X-band using hp-8510 Network Analyzer. For the same composite material using the 
simple mode matching technique the two port S-parameters are calculated as a function 
of complex dielectric constants and thicknesses of each layer. The fminsearch function 
available in the MatLab Optimization Tools is then used to determine complex dielectric 
constants of each layers assuming the thicknesses of each layer are unknown. 

The present approach has been validated using number of composite material 
configurations. The composite material formed by single layers of Garlock, Cork, 
Ceramic, Nomex Felt, Woven Composite B and G, Rubylith, Bakelite, Teflon, Fiber 
Glass, Foam, Nano Material have been tested using the present method. The estimated 
values of dielectric constants and their thicknesses have been found to agree well with the 
true values supplied by the manufacturers. The composite materials formed by placing 
two layers of the above basic materials have been tested using the present method. The 
dielectric constants of individual layers estimated have been found to be in good 
agreement with their values specified by manufacturers. However, for the composite 
material consisting of more than one layer it was assumed that the thicknesses of 
individual layers were known. It has been observed that the present method could not 
estimate the dielectric constants of individual layers correctly when three of more layers 
were used to form a composite material. This is attributed to the fact that the airgap 
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present in the measurement samples could not be accurately modelled in the estimation 

model leading to an incorrect estimation of dielectric constants 
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Appendix A: MatLab code for estimation of S-parameters of multilayer substrate 
placed in a rectangular waveguide. 

% This program calculates S-parameters of 
% N-layer composite material placed in an 
% X-band rectangular waveguide as a function 
% of frequency. Waveguide dimensions (2.29 cm x 1.02 cm) 

% Select 201 frequency points over 8.2-12.4 GHz 
% Frequency band, 
delf = (12.4-8.2)/201; 
forj = 1:201 
fr(j) = 8.2 +(j-l)*delf; 
end 

% Input number of layers 
nl = input('nl = '); 

% Input dielectric constants 
for k = l:2*nl 

% First value is real part of dielectric constant 
% second value is imaginary part of dielectri constant 
x0(k) = input('x0(k) =’); 
end 

% Input location of interfaces 
for k = l:nl+l 

zn(k) = input('zn(k) ='); 
end 

% Set dielectric constants depending upon 
% number of layers. Only 4 layers are 
% included. If nl > 4 then add condition 
% for given nl 


if nl == 1 

er(l) = x0(l) -i*x0(2); 
end 

if nl == 2 

er(l) = x0(l) -i*x0(2); 
er(2) = x0(3) -i*x0(4); 

end 


if nl == 3 

er(l) = x0(l) -i*x0(2); 
er(2) = x0(3) -i*x0(4); 
er(3) = x0(5)-i*x0(6); 

end 
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if nl == 4 

er(l) = xO(l) -i*x0(2); 
er(2) = x0(3) -i*x0(4); 
er(3) = x0(5)-i*x0(6); 
er(4) = x0(7)-i*x0(8); 
end 

nlayer = nl; % number of layers redifined 
n = nlayer+1; % Matrix order 

t =zeros(n,l); % Right handside column matrix initilized 
s =zeros(n); % Admittance matrix initilized 
ak =zeros(n,l); % propagation column matrix 

%Set dielectric constants redifined 

% const is set to zeros (3 by 3) 
const = zeros(n); 

% Dimensions of x-band waveguide in cm 
aa = 2.29; 

bb = 1.02; % guide dimensions 

% Free space dielectric constant 
eO = l.e-9/(36*pi); 

%f = 8.2e9 % frequency of microwave signal 
c = 3.el0; % velocity of light in meters/second 

uO =4*pi*l.e-7 ; % Magnetic permeability of free space 
nO = (uO/eO) A .5; % Free space impedance 
for ifr = 1:201 

freq = fr(ifr); % Frequency in GHz 


akO =2*pi*freq/30; % wave number in free space 
% This loop calculates propagation 
% constants in layers 
for kl = l:nl 

const(kl) = akO*akO*(er(kl))-(pi/aa) A 2; 

if (real(const(kl)) < 0.0) ak(kl) = -i*(-const(kl))\5; 

else ak(kl) = (const(kl)) A .5; 

end 

end 


% akOO and akll are propagation constants 
% for the free space regions of waveguide 
constO = akO A 2-(pi/aa) A 2; 


53 



if (real(constO) < 0.0) akOO = -i*(-constO) A .5; 
else akOO = (constO) A .5; 

end 

akll = akOO; 
for kl = l:nl 

addm(kl) =ak(kl)./(akO*nO); 
end 

addmOO = akOO/(akO*nO); 
addmll = addmOO; 

% Initilize s matrix 
for kl= l:nl+l 
for k2 = l:nl+l 
s(kl,k2) = 0.0; 
end 
end 

% LOOP TO CALCULATE IMPEDANCE MATRIX 
% EXCLUDING FIRST ROW AND LAST ROW 
for kl = 2:nl 

k2 = kl-1; 

argl = ak(k2)*(zn(kl)-zn(kl-l)); 
s(kl.k2) = -addm(k2)/(i*sin(argl)); 
k2 = kl; 

argl = ak(k2-l)*(zn(kl)-zn(kl-l)); 
arg2 = ak(kl)*(zn(kl+l)-zn(kl)); 

s(kl,k2) = addm(k2-l)*cos(argl)/(i*sin(argl)) + addm(kl)*cos(arg2)/(i*sin(arg2)); 
k2 = kl+1; 

arg2 = ak(kl)*(zn(kl+l)-zn(kl)); 
s(kl,k2)= -addm(kl)/(i*sin(arg2»; 
end 

% FIRST ROW OF IMPEDANCE MATRIX 

s(l,l) = addmOO + cos(ak(l)*zn(2))/(i*sin(ak(l)*zn(2)))*addm(l); 
s(l,2) = -addm(l)/(i*sin(ak(l)*zn(2))); 

% LAST ROW OF IMPEDANCE MATRIX 

s(nlayer+l,nlayer) = -addm(nlayer)/(i*sin(ak(nlayer) ... 

*(zn(nlayer+l)-zn(nlayer)))); 
s(nlayer+l,nlayer+l) = addm(nlayer)*cos(ak(nlayer)* ... 
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(zn(nlayer+l)-zn(nlayer)))/(i*sin(ak(nlayer)* 
(zn(nlayer+l)-zn(nlayer)))) + addmll; 


• •• 


% ONLY FIRST ELEMENT OF RIGHT HAND SIDE 
% IS NON ZERO 
for j = l:nlayer+l 
t(j) = 0.0 +i*0.0; 
end 


t(l) =2*addm00; 

% INVERT S AND MULTIPLY BY t 


T = inv(s)*t; 

% FIRST ELEMENT OF T - 1 IS REFLECTION 
% COEFFICIENT AND LAST VALUE OF T is TRANSMISSION 
% COEFFICIENT 

% SINCE MEASURED TRANSIMMSION IS REFERENCED TO 
% INPUT PLANE IT IS MULTIPLIED BY PHASE CORRECTION 


sllc(ifr) =T(1,1)-1; 

s21c(ifr) =T(nl+l,l)*exp(i*akOO*zn(nl+l)); 


end 


% CALCULATING S22 AND S12 


% Inverse the position of interfaces 


if nl == 1 

zn(l) = zn(l); 
zn(2) = zn(2); 
er(l) =x0(l) -i*x0(2); 

end 

if nl == 2 

zn(l) = -zn(3) + zn(3); 
zn(2) = -zn(2) + zn(3); 
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zn(3) = -zn(l) + zn(3); 
er(l) = x0(3) -i*x0(4); 
er(2) = xO(l) - i*x0(2); 

end 

if nl == 3 

zn(l) = -zn(4) + zn(4); 
zn(2) = -zn(3) + zn(4); 
zn(3) = -zn(2) + zn(4); 
zn(4) = -zn(l) + zn(4); 

er(l) = x0(5) -i*x0(6); 
er(2) = x0(4) -i*x0(3); 
er(3) = xO(l) -i*x0(2); 
end 


if n == 4 

zn(l) = -zn(5) + zn(5); 
zn(2) = -zn(4) + zn(5); 
zn(3) = -zn(3) + zn(5); 
zn(2) = -zn(2) + zn(5); 
zn(l) = -zn(l) + zn(5); 
er(l) = x0(7) -i*x0(8); 
er(2) = x0(5) -i*x0(6); 
er(3) = x0(3) -i*xo(4); 
er(4) = xO(l) - i*x0(2); 
end 

nlayer = nl; % number of layers redifined 
n = nlayer+1; % Matrix order 

t =zeros(n,l); % Right handside column matrix initilized 
s =zeros(n); % Admittance matrix initilized 
ak =zeros(n,l); % propagation column matrix 
%Set dielectric constants redifined 

% const is set to zeros (3 by 3) 
const = zeros(n); 
aa = 2.29; 

bb = 1.02; % guide dimensions 
% Free space dielectric constant 
eO = l.e-9/(36*pi); 

%f= 8.2e9 % frequency of microwave signal 
c = 3.el0; % velocity of light in meters/second 

uO =4*pi*l.e-7; % Magnetic permeability of free space 
nO = (uO/eO) A .5; % Free space impedance 
for ifr = 1:201 
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freq = fr(ifr); % Frequency in GHz 
akO =2*pi*freq/30; % wave number in free space 
% This loop calculates propagation 
% constants in layers 
for kl = lsnl 

const(kl) = akO*akO*(er(kl))-(pi/aa) A 2; 

if (real(const(kl)) < 0.0) ak(kl) = -i*(-const(kl)) A .5; 

else ak(kl) = (const(kl)) A .5; 

end 

end 

% akOO and akll are propagation constants 
% for the free space regions of waveguide 
constO = akO A 2-(pi/aa) A 2; 
if (real(constO) < 0.0) akOO = -i*(-constO) A .5; 
else akOO = (constO) A . 5; 

end 

akll = akOO; 

% calculate addmitince for layers: 
for kl = l:nl 

addm(kl) =ak(kl)./(ak0*n0); 
end 

addmOO = ak00/(ak0*n0); 
addmll = addmOO; 

% Initilize s matrix 
for kl= l:nl+l 
for k2 = l:nl+l 
s(kl,k2) = 0.0; 
end 
end 

% LOOP TO CALCULATE IMPEDANCE MATRIX 
% EXCLUDING FIRST ROW AND LAST ROW 
for kl = 2:nl 

k2 = kl-1; 

argl = ak(k2)*(zn(kl)-zn(kl-l)); 
s(kl,k2) = -addm(k2)/(i*sin(argl)); 
k2 = kl; 

argl = ak(k2-l)*(zn(kl)-zn(kl-l)); 
arg2 = ak(kl)*(zn(kl+l)-zn(kl)); 

s(kl.k2) = addm(k2-l)*cos(argl)/(i*sin(argl)) + addm(kl)*cos(arg2)/(i*sin(arg2)); 
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k2 = kl+1; 

arg2 = ak(kl)*(zn(kl+l)-zn(kl)); 
s(kl,k2)= -addm(kl)/(i*sin(arg2)); 


% FIRST ROW OF IMPEDANCE MATRIX 

8(1,1) = addmOO + cos(ak(l)*zn(2))/(i*sin(ak(l)*zn(2)))*addm(l); 
s(l,2) = -addm(l)/(i*sin(ak(l)*zn(2»); 

% LAST ROW OF IMPEDANCE MATRIX 

s(nlayer+l,nlayer) = -addm(nlayer)/(i*sin(ak(nlayer) ... 

*(zn(nlayer+l)-zn(nlayer)))); 
s(nlayer+l,nlayer+l) = addm(nlayer)*cos(ak(nlayer)* ... 
(zn(nlayer+l)-zn(nlayer)))/(i*sin(ak(nlayer)* ... 
(zn(nlayer+l)-zn(nlayer)))) + addmll; 


% ONLY FIRST ELEMENT OF RIGHT HAND SIDE 

% IS NON ZERO 

for j = l:niayer+l 

t(j) = 0.0 -i*0.0; 

end 

t(l) =2*addm00; 

% INVERT S AND MULTIPLY BY t 

T = inv(s)*t; 

% FIRST ELEMENT OF T - 1 IS REFLECTION 
% COEFFICIENT AND LAST VALUE OF T is TRANSMISSION 
% COEFFICIENT 

% SINCE MEASURED TRANSIMMSION IS REFERENCED TO 
% INPUT PLANE IT IS MULTIPLIED BY PHASE CORRECTION 


ref =T(1,1)-1; 

s22c(ifr) = ref*exp(i*ak00*2*zn(nl+l)); 
s!2c(ifr) =T(nl+l,l)*exp(i*akOO*zn(nl+l»; 


end 

ss = load('FT_GRLK.60); 
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delf = (12.4-8.2)/201; 
for j = 1:201 
fr(j) = 8.2 +(j-l)*delf; 
end 

forj = 1:201 

sllr(j) = ss(j,l); 
slli(j) = ss(j,2); 
s21r(j) = ss(j+201,l); 
s21i(j) = ss(j+201,2); 
sl2r(j) = ssO+402,1); 
sl2i(j) = ss(j+402,2); 
s22r(j) = ss(j+603,l); 
s22i(j) = ss(j+603,2); 
end 

plot(fr,sllr,fr.s 1 li,fr,real(sl lc),fr.imag(s 1 lc)) 
ok = input('Is plot ok =’); 
plot(fr,s21r,fr.s21i,fr,real(s21c),fr,iniag(s21c)) 
ok = input('Is plot ok =’); 
piot(fr,s 12r,fr,s 1 2i,fr,real(s 1 2c),f r,imag(s 1 2e )) 
ok = input('Is plot ok ='); 
plot(fr,s22r,fr,s22i.fr,real(s22c),fr,iinag(s22c)) 
ok = input('Is plot ok =’); 


% 
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Appendix B: MatLab code for estimation of dielectric parameters of multilayer 
substrate placed in a rectangular waveguide. 


% This program estimates dielectric parameters of 
% multi-layer substrate . Measured data has to be loaded in 
% meas.dat file 

% X-band rectangular waveguide as a function 
% of frequency. 

%Load measured data file 

ss = load('meas.dat'); 


delf = (12.4-8.2 )/201; 
for j = 1:201 

fr(j) = 8.2 +(j-l)*delf; 

end 

% Input number of layers 
nl = input('nl = '); 

% Input dielectric constants 
for k = l:2*nl 

% First value is real part of dielectric constant 
% second value is imaginary part of dielectri constant 
x0(k) = inputCxO(k) =’); 
end 

% Input location of interfaces 

for k = l:nl+l 

zn(k) = inputCzn(k) ='); 
end 

for j = 1:201 

sllr(j) = ss(j,l); 
sllifj) = ss(j,2); 
s21r(j) = ss(j+201,l); 
s21i(j) = ss(j+201,2); 
sl2r(j) = ss(j+402,l); 
sl2i(j) = ss(j+402,2); 
s22r(j) = ss(j+603,l); 
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s22i(j) = ss(j+603,2); 
end 

% CHECK WHETHER SYSTEM IS LOSSLESS OR LOSSY 

for j = 1:201 
sll = sllr(j)+i*slli(j); 
s21 = s21r(j)+i*s21i(j); 
sl2 = sl2r(j)+i*sl2i(j); 
s22 = s22r(j)+i*s22i(j); 
freq = fr(j); 


xO = [1.0;0.0;1.0;0.0;1.0;0.0]; 

options = optimset('Display', 'iter', 'MaxFunEvals' , 1000, 'TolFun', l.e-08) 
options = optimset('TolX', le-08) 

[x,fval,flag] = fminsearch(@obj ref,xO, options, si I,s21,sl2,s22,freq,nl ... 
,zn); 

erll(j) = x(l); 
erl2(J) = x(2); 
er21(j) = x(3); 
er22(j) = x(4); 
er31(j) = x(5); 
er32(j) = x(6); 

fgl(j) = fval; 
flag 

end 

outmat = [fr.',erll. , ,erl2.']; 
save('tenlpl','outnlat , ,'-ascii , ) 

outmat = [fr.',er21. , ,er22.']; 
save('temp2 , , , outmat', ’-ascii’) 
outmat = [fr.',er31.',er32.']; 
save('temp3', 'outmat', '-ascii') 
outmat =[fr.',fgl.’]; 
save('temp4', 'outmat', '-ascii') 


% 
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% This is object function 

function f = obj_ref( x(),s 1 1 ,s2I ,s 1 2,s22,freq,nl ,zn) 


if nl == 1 

er(l) = xO(l) - i*x0(2); 
end 


if nl == 2 

er(l) = xO(l) -i*x0(2); 
er(2) = x0(3) -i*x0(4); 
end 

if nl ==3 

er(l) = xO(l) -i*x0(2); 
er(2) = x0(3) -i*x0(4); 
er(3) = x0(5) -i*x0(6); 

end 

if nl ==4 

er(l) = xO(l) -i*x0(2); 
er(2) = x0(3) -i*x0(4); 
er(3) = x0(5) -i*x0(6); 
er(4) = x0(7) -i*x0(8); 
end 


nlayer = nl; % number of layers redifined 
n = nlayer+l; % Matrix order 

t =zeros(n,l); % Right handside column matrix initilized 
s =zeros(n); % Admittance matrix initilized 
ak =zeros(n,l); % propagation column matrix 

%Set dielectric constants redifined 

%const is set to zeros (3 by 3) 
const = zeros(n); 
aa = 2.29; 

bb = 1.02; % guide dimensions 
% Free space dielectric constant 
eO = l.e-9/(36*pi); 

%f = 8.2e9 % frequency of microwave signal 
c = 3.el0; % velocity of light in meters/second 

uO =4*pi*l.e-7; % Magnetic permeability of free space 
nO = (uO/eO) A .5; % Free space impedance 
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%for ifr = 1:201 

%freq = fr(ifr); % Frequency in GHz 
akO =2*pi*freq/30; %vvave number in free space 
% This loop calculates propagation 
% constants in three layers 
for kl = l:nl 

const(kl) = akO*akO*(er(kl))-(pi/aa) A 2; 

if (real(const(kl» < 0.0) ak(kl) = -i*(-const(kl))\5; 

else ak(kl) = (const(kl)) A .5; 

end 

end 

constO = akO A 2-(pi/aa) A 2; 

if (real(constO) < 0.0) akOO = -i*(-constO) A .5; 

else akOO = (constO) A . 5; 

end 

akll = akOO; 

% calculate addmitince for layers: 
for kl = l:nl 

addm(kl) =ak(kl)./(ak0*n0); 
end 

addmOO = ak00/(ak0*n0); 
addmll = addmOO; 

% Initilize s matrix 
for kl= l:nl+l 
for k2 = l:nl+l 
s(kl,k2) = 0.0; 
end 
end 

% LOOP TO CALCULATE IMPEDANCE MATRIX 
% EXCLUDING FIRST ROW AND LAST ROW 
for kl = 2:nl 

k2 = kl-1; 

argl = ak(k2)*(zn(kl)-zn(kl-l)); 
s(kl,k2) = -addm(k2)/(i*sin(argl)); 
k2 = kl; 

argl = ak(k2-l)*(zn(kl)-zn(kl-l)); 
arg2 = ak(kl)*(zn(kl+l)-zn(kl)); 

s(kl,k2) = addm(k2-l)*cos(argl)/(i*sin(argl)) + addm(kl)*cos(arg2)/(i*sin(arg2)); 
k2 = kl+1; 
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arg2 = ak(kl)*(zn(kl+l)-zn(kl)); 
s(kl,k2)= -addm(kl)/(i*sin(arg2)) ; 
end 


% FIRST ROW OF IMPEDANCE MATRIX 

s(I,l) = addmOO + cos(ak(l)*zn(2))/(i*sin(ak(l)*zn(2)))*addm(l); 
s(l,2) = -addm(l)/(i*sin(ak(l)*zn(2)»; 

% LAST ROW OF IMPEDANCE MATRIX 

s(nlayer+l,nlaver) = -addm(nlayer)/(i*sin(ak(nlayer) ... 

*(zn(nlayer+l)-zn(nlayer)))); 
s(nlayer+l,nlaver+l) = addm(nlayer)*cos(ak(nlayer)* ... 
(zn(nlayer+l)-zn(nlayer)))/(i*sin(ak(nlayer)* ... 
(zn(nlayer+l)-zn(nlayer)))) + addmll; 


% ONLY FIRST ELEMENT OF RIGHT HAND SIDE 
% IS NON ZERO 
for j = l:nlayer+l 
t(j) = 0.0+i*0.0; 
end 

t(l) =2* addmOO; 

% INVERT S AND MULTIPLY BY t 


T = inv(s)*t; 


% FIRST ELEMENT OF T - 1 IS REFLECTION 
% COEFFICIENT AND LAST VALUE OF T is TRANSMISSION 
% COEFFICIENT 

% SINCE MEASURED TRANSIMMSION IS REFERENCED TO 
% INPUT PLANE IT IS MULTIPLIED BY PHASE CORRECTION 


sllc =T(1,1)-1; 

s21c =T(nl+l,l)*exp(i*akOO*zn(nl+l)); 

f = real(sllc-sll)*real(sllc-sll) + imag(sllc-sll)*imag(sllc-sll) . 
+real(s21c-s21)*real(s21c-s21) + imag(s21c-s21)*imag(s21c-s21); 


% CALCULATING S22 AND S12 
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if nl ==1 

znl(l) = zn(2) - zn(2); 
znl(2) = zn(2) - zn(l); 
end 

if nl == 2 

znl(l) = zn(3) - zn(3); 
znl(2) = zn(3) - zn(2); 
znl(3) = zn(3) - zn(l); 
end 


if nl == 3 

znl(l) = zn(4) - zn(4); 
znl(2) = zn(4) - zn(3); 
znl(3) = zn(4) - zn(2); 
znl(4) = zn(4) - zn(l); 
end 


if nl == 4 

znl(l) = zn(5) - zn(5); 
znl(2) = zn(5) - zn(4); 
znl(3) = zn(5) - zn(3); 
znl(4) = zn(5) - zn(2); 
znl(5) = zn(5) - zn(l); 
end 


if nl == 1 

er(l) =xO(l) - i*x0(2); 
end 

if nl == 2 

er(2) = xO(l) - i*x0(2); 
er(l) = x0(3) - i*x0(4); 

end 

if nl ==3 

er(l) = x0(5) -i*x0(6); 
er(2) = x0(3) -i*x0(4); 
er(3) = xO(l) - i*x0(2); 
end 

if nl ==4 

er(l) = x0(7) -i*x0(8); 
er(2) = x0(5) -i*x0(6); 
er(3) = x0(3) -i*x0(4); 
er(4) = xO(l) - i*x0(2); 

end 


nlayer = nl; % number of layers redifined 
n = nlayer+1; % Matrix order 

t =zeros(n,l); % Right handside column matrix initilized 
s =zeros(n); % Admittance matrix initilized 
ak =zeros(n,l); % propagation column matrix 
%Set dielectric constants redifined 

%const is set to zeros (3 by 3) 
const = zeros(n); 
aa = 2.29; 

bb = 1.02; % guide dimensions 
% Free space dielectric constant 
eO = l.e-9/(36*pi); 

%f = 8.2e9 % frequency of microwave signal 
c = 3.el0; % velocity of light in meters/second 

uO =4*pi*l.e-7; % Magnetic permeability of free space 
nO = (uO/eO) A .5; % Free space impedance 
%for ifr = 1:201 

% freq = fr(ifr); % Frequency in GHz 
akO =2*pi*freq/30; % wave number in free space 
% This loop calculates propagation 
% constants in three layers 
for kl = l:nl 

const(kl) = akO*akO*(er(kl))-(pi/aa) A 2; 

if (real(const(kl)) < 0.0) ak(kl) = -i*(-const(kl)) A .5; 

else ak(kl) = (const(kl)) A .5; 

end 

end 

%ak 

%ok = input('ok = ') 

% akOO and akll are propagation constants 
% for the free space regions of waveguide 
constO = akO A 2-(pi/aa) A 2; 
if (real(constO) < 0.0) akOO = -i*(-constO) A .5; 
else akOO = (constO) A .5; 

end 

akll = akOO; 

%ak0 

% calculate addmitince for layers: 
for kl = l:nl 

addm(kl) =ak(kl)./(ak0*n0); 
end 
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%addm 

%ok = inputCok =’) 
addmOO = ak00/(ak0*n0); 
addmll = addmOO; 

% addmOO 
% addmll 
%ok = input('ok =') 

% Initilize s matrix 
for kl= l:nl+l 
for k2 = l:nl+l 
s(kl,k2) = 0.0; 
end 
end 

% LOOP TO CALCULATE IMPEDANCE MATRIX 
% EXCLUDING FIRST ROW AND LAST ROW 
for kl = 2:nl 

k2 = kl-1; 

argl = ak(k2)*(znl(kl)-znl(kl-l»; 
s(kl,k2) = -addm(k2)/(i*sin(argl)); 
k2 = kl; 

argl = ak(k2-l)*(znl(kl)-znl(kl-l)); 
arg2 = ak(kl)*(znl(kl+l)-znl(kl)); 

s(kl,k2) = addm(k2-l)*cos(argl)/(i*sin(argl)) + addm(kl)*cos(arg2)/(i*sin(arg2)); 
k2 = kl+1; 

arg2 = ak(kl)*(znl(kl+l)-znl(kl»; 
s(kl,k2)= -addm(kl)/(i*sin(arg2)); 
end 

% FIRST ROW OF IMPEDANCE MATRIX 

s(l,l) = addmOO + cos(ak(l)*znl(2))/(i*sin(ak(l)*znl(2)))*addnt(l); 
s(l,2) = -addm(l)/(i*sin(ak(l)*znl(2))); 

% LAST ROW OF IMPEDANCE MATRIX 

s(nlayer+l,nlayer) = -addm(nlayer)/(i*sin(ak(nlayer) ... 

*(znl(nlayer+l)-znl(nlayer)»); 
s(nlayer+l,nlayer+l) = addm(nlayer)*cos(ak(nlayer)* ... 
(znl(nlayer+l)-znl(nlayer)))/(i*sin(ak(nlayer)* ... 
(znl(nlayer+l)-znl(nlayer)))) + addmll; 

%s 

%ok = inputCok =') 

% form cpmplex matrix 
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% ONLY FIRST ELEMENT OF RIGHT HAND SIDE 
% IS NON ZERO 
for j = l:nlayer+l 
t(j) = 0.0 + i*0.0; 
end 

t(l) =2*addm00; 

% INVERT S AND MULTIPLY BY t 


T = inv(s)*t; 


% FIRST ELEMENT OF T - 1 IS REFLECTION 
% COEFFICIENT AND LAST VALUE OF T is TRANSMISSION 
% COEFFICIENT 

% SINCE MEASURED TRANSIMMSION IS REFERENCED TO 
% INPUT PLANE IT IS MULTIPLIED BY PHASE CORRECTION 


ref =T(1,1)-1; 

s22c = ref*exp(i*ak00*2*znl(nl+l»; 
sl2c =T(nl+l,l)*exp(i*akOO*znl(nl+l)); 

f = f+ real(s22c-s22)*real(s22c-s22) +imag(s22c-s22)*imag(s22c-s22) 
+ real(s!2c-sl2)*real(sl2c-sl2) + imag(sl2c-sl2)*imag(sl2c-sl2); 
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